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ABSTRACT 

The high efficiency of converting kinetic energy into gamma-rays estimated 
with late-time afterglows in Gamma-Ray Burst (GRB) phenomenon challenges 
the commonly accepted internal-shock model. However, the efficiency is still 
highly uncertain because it is sensitive to many effects. In this Letter we study 
the sideways expansion effect of jets on estimating the efficiency. We find that 
this effect is considerable, reducing the efficiency by a factor of ~ 0.5 for typical 
parameters, when the afterglow data ~ 10 hr after the GRB trigger are used 
to derive the kinetic energy. For a more dense circumburst medium, this effect 
is more significant. As samples, taking this effect into account, we specifically 
calculate the efficiency of two bursts whose parameters were well constrained. 
Almost the same results are derived. This suggests that the sideways expansion 
effect should be considered when the GRB efficiency is estimated with the late 
afterglow data. 

Subject headings: gamma ray: bursts — gamma ray: observations 

1. Introduction 

The efficiency [r]^) of converting kinetic energy into 7-rays (E^) is one of the most 
important physical parameters of gamma-ray burst (GRB) phenomenon (Zhang et al. 2007, 
hereafter Z07). In the conventional internal+external shock GRB model (Rees & Meszaros 
1994; Kobayashi, Piran & Sari 1997; see recent review by Zhang & Meszaros 2004), 77^ is 
~ 1% — 5% (e.g., Kobayashi, Piran & Sari 1997; Kumar 1999), and even under some extreme 
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assumptions, such as the extremely inhomogeneous velocity of the ejecta shells, r)j is only 
40% (Kobayashi, Piran & Sari 1997, Kobayashi & Sari 2001). However, the r]^ estimated 
with the late afterglow (say, 10 hr after the GRB trigger) are much higher (> 50%) than the 
model prediction (Lloyd- Ronning & Zhang 2004, hereafter LZ04). By considering the inverse 
Compton effect. Fan & Piran (2006, hereafter FP06) derived r]^ = 1% ~ 89%. More recently, 
Z07 revisited this issue for 32 GRBs detected by Swift with the early X-ray afterglow data. 
They derived the kinetic energy (isotropic-cquivalent) at the end of the shallow decay phase 
(tf)) and at the deceleration time of the fireball (tdec) , ^-nd they found that r)^ estimated with 
the kinetic energy at tb is < 10% for most of the GRBs and that the rjj estimated with 
the kinetic energy at tdec varies from a few percent to > 90%. These results challenge the 
conventional internal-shock model. 

The difficulty of measuring rj^ is the derivation of the kinetic energy (Ek) of the af- 
terglows, which sensitively depends on many physical parameters of the burst itself and its 
environment (Z07). Measurement of the Ek with early afterglow data seems more reliable 
(Z07). However, it is difficult to determine the onset of afterglow from the observation. 
Furthermore, most of the X-ray early afterglows detected by Swift show complex features 
with irregular flares (Zhang et al. 2006; Nousek et al. 2006; Panaitescu et al. 2006; Granot 
et al. 2006). Estimating the kinetic energy with them can also introduce some uncertain- 
ties. Zhang et al. (2007) systematically analyzed the prompt and early afterglow emission 
of 32 Swift bursts and found that the GRB efficiency derived with early afterglow is highly 
variable. This result — if not intrinsic — seems to prove the above viewpoint. The derivation 
of Ek with the late afterglow data, on the other hand, involves some uncertainties due to 
the energy loss, energy injection, and jet sideways expansion effect. The advantage of the 
late afterglow is that this phase is steady, generally showing a smooth power-law decay. 
By taking the energy loss, the injection, and the jet expansion effect into account, the late 
afterglow may also be used for estimating r]^. 

It is now believed that GRBs arise from jets collimated into a small angle. To derive the 
actual radiative and kinetic energy from the observed fluence, the collimation correction has 
to be considered. In previous works, the isotropic- equivalent E^ and E^ were used to derive 
the efficiency without considering this geometrical correction, which is only valid when the 
effect of the jet expanding sideways is ignored. However, the GRB jet evolves and expands 
sideways with time. Although this effect is not notable at a very early stage (Huang et al. 
2000), it may affect the estimation of the efficiency 77^, when one uses late afterglows at ~ 10 
hr to calculate the efficiency. In this Letter, we investigate the jet sideways expansion effect 
on estimating the GRB efficiency. Our theoretical analysis is present in §2. Cases studies 
for some typical GRBs are present in §3. Discussion and conclusion are presented in §4. 
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2. Theoretical Analysis 



The 1]^ is defined as 



'^'^ (1) 

Assuming fb^^ and fb,k are the beaming factors during the prompt and the afterglow phases, 
respectively, the corrected efficiency r]^ is given by 

E,h,k + E,h,, / + (!-/)//,' 

where / = fb,k/ fb,-y It can be seen that when jet angle evolution is considered, it is equivalent 
to increasing the E^ by a factor of /. For a double-sided jet, fb^-y = 1 — cos{9^), and 
fb,k = 1 — cos{9k), where 9^ and 9k are the opening angles during the prompt and the 
afterglow phases, respectively. 

The hydrodynamic evolution of the jet has been studied by a number of authors (e.g., 
Rhoads 1997, 1999; Panaitescu & Meszaros 1999; Sari, Piran, & Halpern 1999; Moderski, 
Sikora, & Bulik 2000; Huang et al. 2000; Kumar & Panaitescu 2000). We adopt the evolution 
equation group presented in Huang et al. (2000), which describes the overall evolution of 
relativistic GRB ejecta from the an ultra-relativistic phase to a non-relativistic phase. The 
equations are quoted as follows: 

^ = /3c7(7+v/y^), (3) 
^ = 27ri?2/5c7(7 + ^7^ - 1)(1 - cos9)nmp, (4) 



d9 _ c,(7+ V^^) 
dt~ R ' 

d-f _ _27ri?2/5c7(7 + ^7^ - 1)(72 - 1)(1 - cos9)nmp 

dt ~ Mej + em + 2(1 - e)-fm ' ^ ^ 

where R is the radial coordinate in the burst source frame, 7 is the bulk Lorentz factor 
of the ejecta [/3 = (7^ - lY^^j], m is the swept-up mass, n is the number density of the 
surrounding interstellar medium (ISM), rrtp is the mass of the proton, 9 is the half-opening 
angle of the jet, Cs = {[(7(7 — 1)(7 — 1))]^^^/[1 +7(7 ~ is the comoving sound speed, 

7 = (47 + l)/(37) is the adiabatic index, and e is the radiation efficiencyll] Mej is the ejecta 
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mass that is defined by Eiso{l — cosOq) — ^yoMejc^, where Eiso, Oo, and 70 are the burst energy 
(isotropic- equivalent), the initial half-opening angle of jet, and the initial Lorentz factor, 
respectively. 

We first adopt typical parameters as follows in order to study the jet sideways effect: 
70=300, Eiso— 10^^ ergs, 9o—0.1 radians, and n — 1 cm~^. In our calculation, we take 9j to 
be the initial half-opening angle of the jet and 9k as the one at 10 hr. With these typical 
parameters, we derive the evolution of the jet opening angle, which is shown in Figure 1. We 
find that 9^ ~ 0.14 at ~ 10 hr. Therefore, / ~ 2, and hence = ri^/{2 — ri^). This indicates 
that the lower the efficiency, the more significant the reduction is on the efficiency from this 
effect. For a burst with r]^ = 50%, 77^ = 2ri^/3, while for rj^ = 10%, rj^ = r)^/2. The initial 
parameters, E^so and n, vary between bursts. Figure 1 also shows the evolution of the jet 
opening angle for various values of Eiso a-^d n; rf^ as a function of the initial parameters is 
shown in Figure 2. It is found that the correction to 77-y tends to be more significant for a 
burst with smaller ^0, lower Eiso-, and a denser ambient medium. 



3. Case Studies 

We take two pr(^Swift bursts, GRB 980703 and GRB 000926, as examples. The micro- 
physics parameters of the two bursts have been well constrained with high-quality broadband 
afterglow data (Yost et al 2003). For calculating a more accurate hydrodynamic evolution, 
the radiative loss during the early afterglow phase cannot be neglected. Since the energy of 
the accelerated electrons behind the blast wave is lost both through the synchrotron radia- 
tion (inverse Compton cooling is neglected) and through the expansion of the fireball, the 
radiative efficiency is the ratio of the synchrotron power to the total power, namely (Dai & 
Lu 1998,1999b) 

^-^ejjf^„ (7) 

syn ' ex 

where eg is the fraction of shock energy given to the electrons, t'sy^ = OTrmgc/ {crTB'^je,min) 
is the synchrotron cooling time, with ax the Thompson cross section and 7e,mm — ^e{l ~ 
l)mp{p — 2)/me{p — 1) -|- 1 the minimum Lorentz factor of the random motion of electrons 
in the comoving frame, and t'^^ — R/ (7c) is the comoving frame expansion time. With the 
numerical factor derived by FP06, we derive 

Ek = lO^^ergs ^(p) Lf,f (i±^)(2--)/(-+2)6^(!-^)/(-+^)4(V)/(-+^\ (8) 
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Fig. 1. — Evolution of the jet half-opening angle. The vertical solid line marks the position 
of t=10 hr. The thick solid line is for the typical parameters. The dotted lines are for the 
typical parameters with n = 10,000, 100, and 0.01 cm ~^ from the top to the bottom. The 
dashed lines are for the typical parameters with Eigo = 10^^, 10^^, and 10^^ ergs from the 
top to the bottom. Note that the hnes for n — 100 cm~^ and Eiso — 10^^ ergs are almost 
superposed. 
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Fig. 2. — The rj^ in some extent of the initial parameters, which is normalized to 77 = 0.1 
and the initial parameters to the typical values. When one parameter varies, the others are 
taken as the typical values. 
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where 



A{p) = { 



1.35 X 10^ X (1 -p/2) 




0.42(2-3p)/4 6 X i0ii)(p-i)/2 C^-^^ [(2.42 x iQ^syi-pm _ (4 34 x ioi6)(i-p/2)] 



Here p is the spectral index of the electron distribution, and Cp = 13(p — 2)/3(j9 — 1), 
Ri ~ [t(10/i)/T9o]^^'^'=/^^; is a factor that accounts for radiative losses during the first 10 
hr following the prompt phase (Sari 1997, LZ04), Lx,46 = Lx/lO"^^ is the isotropic X-ray 
afterglow luminosity at 10 hr (using the value given by Berger et al. 2003), and is the 
fraction of shock energy given to the magnetic field (the inverse Compton effect is neglected). 
Ek is the isotropic kinetic energy at i = 10 hr. 

Assuming the jet opening angle at 10 hour is ^^lo, the initial kinetic energy of cjccta is 
Eq = Ek{l — cosOiq). Wc thus obtain a boundary condition that when t = tio, 9 = Oiq. 
In addition, in the standard model, there is a light-curve break occurring at the time (tj) 
when the bulk Lorentz factor of the shock has slowed to 7 ~ (Rhoads 1997; Sari et al. 
1999). This can be taken as another boundary condition; namely, when t — tj, ^ — 6~^. 
Note that the second condition may depend on the assumed initial Lorentz factor. However, 
we find that the calculated evolutions of the Lorentz factor with different assumed initial 
values (from 100 to 500) arc almost the same during the late time (>1000 s with the typical 
parameters), which is understandable. At the late time, the mass of the jet is dominated by 
the swept-up mass, so the initial mass of ejecta, which is defined by the initial Lorentz factor 
with a given initial kinetic energy, weakly affects the evolution of the Lorentz factor of a jet. 
We adjust 6iq and 6q in a proper range until the above two boundary conditions are satisfied. 
Then the initial jet angle and the initial kinetic energy are derived, and the efficiency can 
be calculated. Our results are hsted in Table 1. It can be seen that the corrected efficiency 
is significantly lower than that derived by LZ04 and FP06. Note that the E^ derived in this 
analysis is different from LZ04 and FP06 due to the different microphysics parameters used. 
If we use the kinetic energy that we obtain, the efficiencies of GRB 980703 and GRB 000926 
derived with equation (1) are 0.135 and 0.103, respectively, while with equation (2), they 
are only 0.079 and 0.052, respectively. 



We have investigated the correction on the estimate of 77^ by considering the jet sideways 
effect. Our numerical study shows that for a stronger burst, a burst with a narrower jet 
opening angle, or a burst in a more dense environment (Dai & Lu 1999a), this effect tends to 
be more significant. Using the typical parameters, the corrected efficiency may be reduced 



4. 



Discussions and conclusion 
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by a factor of ~ 0.5 for a GRB with typical parameters. We use two pie-Swift bursts, GRB 
980703 and GRB 000926, as detailed case studies, further confirming our numerical analysis 
results. 

Several authors argued that the efficiencies of some Swift bursts with long time X-ray 
flattening are as high as ~75%- 90% (Zhang et al. 2006; Nousek et al. 2006; loka et al. 
2006; Granot et al. 2006) and that the internal-shock model for GRBs meets the challenge 
of a so-called efficiency crisis (e.g., loka et al. 2006). As mentioned in §1, the estimate 
of the GRB efficiency is strongly affected by many effects and depends on unobservable 
physical parameters. As we show in this analysis, the small jet sideways effect may result 
in a considerable reduction of the efficiency. More recently, Zhang et al. (2007) estimated 
the kinetic energy with the data of early afterglows, such as the deceleration time (when the 
blast wave is decelerated). From the end of GRB to the deceleration time, the jet sideways 
expansion is insignificant, and thus this effect on the estimate of the GRB efficiency can be 
negligible. However, they found that the derived GRB efficiency is highly variable from a 
few percent to >90% (at the deceleration time). If these results are intrinsic but not caused 
by a parameter selection effect or by the uncertainties of early afterglows, they may refresh 
our understanding of the GRB phenomenon, and they may require an improvement in the 
conventional shock model. 

We thank Y. Z. Fan for helpful discussion and suggestions. We also acknowledge the 
referee for constructive comments and suggestions. This work is supported by the National 
Natural Science Foundation of China (grants 10443003 and 10573030) and the Natural Sci- 
ence Foundation of Yunnan (2003 A0025Q). 
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Table 1: The parameters and the calculated GRB efficiency. is the isotropic radiative 
energy from Bloom et al. (2003), microphysics parameters, n and tj from Yost et al. (2003). 
Here we assume 70 = 300. a and b are the efficiency from LZ04 and FP06, respectively. 



GRB 


Ey 

(lO^^ergs) 


^0 




es 


P 


(lO^^ergs) 


n 

(cm~^) 


(days) 


^7 
(a,b) 


980703 


6.01 


8.37° 


0.27 


0.0018 


2.54 


3.85 


28 


3.4 


0.08 
(0.71,0.21) 


000926 


27.97 


4.59° 


0.15 


0.022 


2.79 


24.3 


16 


2.6 


0.05 
(0.74,0.23) 



